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We study the influence of quantum density fluctuations in ultracold atoms in an optical lattice 
on the scattering of matter waves. Such fluctuations are characteristic of the superfluid phase 
and vanish due to increased interactions in the Mott insulating phase. We employ an analytical 
treatment of the scattering and demonstrate that the fluctuations lead to incoherent processes, 
which we propose to observe via decoherence of the fringes in a Mach-Zender interferometer. In 
this way we extract the purely coherent part of the scattering. Further, we show that the quantum 
density fluctuations can also be observed directly in the differential angular scattering cross section 
| for an atomic beam scattered from the atoms in a lattice. Here we find an explicit dependence of 

the scale of the inelastic scattering on the quantum density fluctuations. 
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I. INTRODUCTION 



The atoms in a weakly interacting Bose Einstein condensate flow as a superfluid over the crystalline potential 
landscape formed by a periodic pattern of light intensity. Individual atoms in this system can be spatially delocalized 
over the breadth of the optical lattice, and, as a consequence, fluctuations in the number of atoms at individual lattice 
sites occur even at zero temperature. As the repulsion between the atoms grows, they are driven into a Mott insulating 
phase, in which they maximally isolate themselves and localize to a single lattice site [1-3]. This rigid crystal, which 
reflects the behavior of tightly bound electrons hopping in a solid, develops as the atoms spontaneously lock to fixed, 
Q integer occupation of each lattice site, and the fluctuations vanish. 

Our objective is to study the influence of these quantum mechanical fluctuations on the scattering of matter 
waves, and to illustrate the susceptibility of the scattering to this purely quantum mechanical effect, as well as the 
Ch ■ consequent impact of the transition from non-interacting single particle physics to strongly interacting, many body 
physics. Ultracold atoms in an optical lattice are an ideal context in which to examine this behavior because they allow 
i^i . controlled investigation of the transition between single and many body physics, due to the experimental tunability 
of the physical parameters of the crystal: the tunneling rate between lattice sites, the interaction strength between 
, atoms in the lattice and the geometry of the lattice [3-5] . 

It has been proposed that light scattering from atoms in an optical lattice, situated within an optical cavity, will 
be sensitive to the on-site number statistics [6]. Moreover, recent experiments with high resolution optical systems 
have demonstrated images resolving features on the scale of a single lattice site both for atoms collected from a 
magneto-optic trap [7, 8] and for atoms in a Bose- Einstein condensate [5, 9]. Due to the interactions of the atoms 
with the light, however, only the parity of the number of atoms at each site can be directly observed; nonetheless, the 
fluctuations in the parity display a dependence on the many body phase of the atoms in the lattice. 

The scattering cross section of ultracold atoms in an optical lattice seen by a matter wave probe is sensitive to the 
quantum many body phase in the lattice [10]. We seek to understand the specific influence of the exotic properties 
of these materials on the images produced by matter waves and to employ this knowledge to make non-destructive 
measurements. Here, we study the effect of the purely quantum mechanical density fluctuations at individual lattice 
sites and show that these zero-temperature fluctuations can be directly probed using a matter wave. 

Fluctuations in the density distribution seen by an imaging matter wave are expected to produce a dephasing effect 
and lead to incoherent scattering. By examining the overlap of the matter wave scattered from the atoms in an 
optical lattice with a coherent reference beam, in an interferometric configuration, we will be able to explicitly see 
the decohering effect of quantum density fluctuations on the interference pattern, and to determine the size of the 
density fluctuations. Subsequently, we will show that this incoherent scattering can even be observed in the angular 
differential scattering cross section, and that the scale of the fluctuations can be obtained directly by examining the 
inelastic part of the scattering in an appropriate regime of probe energy and lattice depth. 

The matter wave probe is a free particle with mass m, initial wave vector feo and Hamiltonian Hp = p 2 /2m, which 
does not interact with the lattice light. The scattering target is well modeled by the Bose Hubbard Hamiltonian, 
describing N atoms confined to the lowest band of a lattice with Nr, sites, 

H B a = -J ^2 a R aR> + \uy^n R (n R - 1). (1) 

(R,R') R 
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FIG. 1. Scattering configuration in which the probe matter wave is incident from the left upon a one-dimensional lattice. A 
detector in the far-field measures the flux of probe atoms scattered in the direction 6. 



For a sufficiently large lattice, in which the effects of the edges of the lattice are negligible, we may impose periodic 
boundary conditions on both the atoms in the lattice and the transverse dimensions of the probe. We may freely 
choose our coordinate system so that the z-axis is aligned with k . The initial probe wave function is then given by 

p ik z 

<"*■> " 7531' (2) 

where fco is the magnitude of feo, and L is the length of the lattice in the transverse dimensions. 

As we wish our probe to weakly interact with the atoms in the lattice, and to avoid interband excitations, we will 
insist on low energy probes, for which s-wave scattering from the atoms in the target is dominant, and we may treat 
the interaction between the probe and the atoms in the lattice as a pseudopotential with scattering length, a s . The 
scattering interaction is then given by [11], 

V = J2v 5(r-rj), (3) 

3 

where the operators r and rj give the positions of the probe and the j th lattice atom, respectively, and Vq = 2Trh 2 a s /m. 
The full Hamiltonian for the scattering interaction is H = Hp + Hbh + V. 



II. INTERFERENCE MEASUREMENT OF DENSITY FLUCTUATIONS 



A standard scattering configuration, as depicted in Fig. 1, in which the angular differential scattering cross section 
is measured, does not specifically distinguish coherent from incoherent scattering. The possibility that disturbances 
are created in the target, such as long wavelength phonons, exists despite the presence of evidence in the scattering 
pattern for coherent scattering processes, such as Bragg peaks. Insofar as these disturbances do not completely localize 
the probe, interference of its wave function scattered from separated points in the lattice can persist. 

The situation is starkly different if the probe is coherently split into two beams before scattering, so that there is 
a reference beam that does not interact with the target. Interference between the atom scattered from the target 
and the reference beam will then take place only if the probe scatters without disturbing the target at all. We will 
confirm the incoherent effect of the quantum density fluctuations in the lattice on the scattered probe specifically by 
employing the Mach-Zender configuration shown in Fig. 2 [12, 13]. 

In this arrangement, the upper arm scatters from the sample of cold atoms in the optical lattice and then interferes 
downstream with the lower arm. The interference pattern that forms there will be due solely to the coherent part of 
the scattering, in which the state of the atoms in the lattice is undisturbed [14, 15]. In practice, the narrow angular 
acceptance of the region of overlap, where the detector is placed, also restricts detection to the essentially forward 
scattered atoms [16, 17]. The forward scattered beam may, however, acquire a phase shift and be attenuated due 
to incoherent scattering. Both of these quantities, the phase shift and attenuation, can be directly observed in the 
interference pattern that forms. Here, we are particularly interested in analyzing the attenuation, as this indicates the 
amount of incoherent scattering, and, as we demonstrate below, depends on the zero-temperature density fluctuations. 

In order to determine this phase shift and attenuation, we consider the S-matrix for the multi-particle system, 
which connects the initial many body state of the probe and target with the asymptotic output state, \ip), long after 
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FIG. 2. Schematic of the Mach-Zender interferometer arrangement in which the probe beam is incident on a diffraction grating 
from the left. The th and 1 st order diffracted beams form the lower and upper arms of the interferometer, respectively. The 
second diffraction grating causes the arms of the interferometer to overlap downstream at the position of the screen. A target 
is placed in the path of the upper arm only, so that only the coherent scattering will interfere with the lower, reference arm. 



the scattering interaction has ceased [18], 

|V> = £|fc ,no>. (4) 

The initial wave vector of the probe is labeled fcrj, arid |no) is the many body ground state of the N atoms in the 
lattice. The operator S can be expanded generically in terms of the total energy E of the N + l atoms, the transition 
operator T and the Hamiltonian H , which excludes interactions [19], 

S= l-2mS(E-H ) T. (5) 

For the system we are discussing, the non- interacting Hamiltonian is given by Hq = Hp + Hbu, and the transition 
operator is given by [20], 

f = V + lim V 1 ; f. (6) 

*->o E Q + E na -H P - F BH + it 



The operator 



(Eq + E no -Hp- H BH + ie) 1 (7) 



is the free Green's operator for the many body system. Eq is the initial energy of the probe, and E no is the energy 
of the target ground state. The state of the N+l particle system that results from the scattering, in terms of these 
operators, is 

\ij>) = \k ,n ) - 2m5(E - #„) f \k , n ) . (8) 

The first term is the unscattered initial state. In the Mach-Zender configuration, this term is coherent with the 
reference arm. The second term involves both coherent and incoherent contributions. We can facilitate extraction 
of the coherent contribution by expanding in terms of eigenstates of the non-interacting probe and target. This is 
accomplished by interposing a complete set of probe and target states between the delta function and the operator T. 
We may then exclude contributions to the scattering that create disturbances in the target or deflect the probe out 
of the forward direction. Choosing the direction of k to align with the z-axis of our coordinate system and recalling 
that the probe is subject to periodic boundary conditions in the transverse dimensions, the expansion yields 

k z ,n\T\k ,n a ) . (9) 

kx 1 ky ¥1 

Any scattering interaction which changes the state of the target, such that n ^ no, would cause the probe to decohere 
completely, so that we may exclude all such terms from the sum over target states. As we mentioned above, the 
small angular acceptance of the detector requires that we also exclude deflection of the probe in the x — y plane. The 
remaining integral over k z may then be performed immediately. The coherent part of the system after scattering is 
given by 

(27T77J - \ 

1 - (ka,n \T\ko,n ) j \k ,n ) , (10) 
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where, for simplicity, we write |fco) as short hand for |0,0, fco), in which the components of the probe wave vector in 
the x and y directions are zero. Eq. (10) does not yet incorporate the specific details of our target, and is written 
here in general for any scattering matrix T and initial target state \no). The diagonal T-matrix element of the initial 
state is the crucial quantity that determines both the phase shift and the contrast of the interference fringes. We may 
identify the complex amplitude of the coherent wave, given in Eq. (10), with the expansion of an exponential when 
the strength of the interaction Vb between the probe and the atoms in the target is weak [21]. The coherent amplitude 

can then be separated into a part which corresponds to a phase shift, <j> — — |f^Rc (jf)^ , and an attenuation, 

£ = — "^S^Im (jT)^, where both expectation values are taken in the initial state, \k ,n ). The coherent part of the 
scattering in terms of these quantities is 

IV'coh) = (1 + %tf> - |*o, n ) « e'*-« \k ,n ) . (11) 

We may read off from this expression that the interference pattern that forms will be shifted by <j> and the fringes will 
have a contrast e~"> [14, 15]. 

It is appropriate to consider the case of weak interactions as we seek to probe the target without disturbing its 
state. The values of <ft and £ will then be dominated by the lowest-order terms in a Born expansion of the operator 
T, defined by Eq. (6). To first order in V , the transition operator coincides with the interaction potential (T m V), 
so that the T-matrix element becomes (ho, no \ V \ ko, n ). This quantity is purely real, and gives us the leading order 
phase shift of the interference pattern. To first order in Vb, however, there is no attenuation, and the scattering is 
purely coherent. Loss of contrast in this system arises in the second order correction, T w V + VG Q N+1 ^V. The second 

order contribution to the T-matrix element, (fc , n \ VG Q N+1>> V \k , no), is in general complex-valued. The real part of 
it will contribute a small correction to the phase shift (j), and the imaginary part gives the leading order approximation 
of£. 

In order to facilitate our intuition of the relationship between <fi and £ and the moments of the density distribution 
of the atoms in the lattice, we will express the interaction potential in second quantized form, in which it is given by 

V = V n{r ). (12) 

Here n(r) — xft (r)ip(r) is the density of the target at the position r. Inasmuch as the interaction potential depends 
linearly on n(r), an expansion of the T-matrix in a Born series, in orders of Vb, will contain terms that depend upon 
a corresponding moment of the density distribution. For this reason, 0, whose leading contribution is linear in Vb, 
will depend solely on the average density of the target. Likewise, £, which is second order in Vb, will depend on the 
second moment of the density distribution of atoms in the target, and therefore on the fluctuations. 

Let us now explicitly evaluate the phase shift and decoherence of the probe. As we explained above, the phase shift 
requires that we evaluate the diagonal T-matrix element in the first Born approximation, which yields 

V f 

(fc ,no|T|fco,n ) ~ l^Jl d 3 r (n \ h(r) \n ) . (13) 

We find that, for any |no), to first order in the interaction strength, in which only single scattering events are included 
in the expansion of the T-matrix, the coherently scattered wave depends only on the average density of the many 
body target, 

l^coh) = exp (^-i —j^^ J l fc o,«o) • (14) 

The approximation of the phase shift given by Eq. (14) is proportional to the column density, p = N/L 2 , of the atoms 
in the target. This result is not restricted to the case of cold atoms in an optical lattice, and is known to determine, 
for example, the phase shift due to scattering from a thermal gas [15, 16, 22, 23]. 

By allowing terms that are second order in Vb, we enable scattering channels in which the probe scatters from two 
different atoms within the target, and thereby becomes sensitive to density fluctuations. The coherent part of the 
scattered wave to second order is, 

IV'coh) = (^l-^^( K (ko,no\V\ko 7 no) + (ko,no\VG { N+1) V\ko,no)^ \k ,n ) . (15) 

As we argued above, the imaginary part of the second order term gives us the exponent in the fringe contrast to 
leading order, so that we must evaluate 

2nm 



(k ,n \VGi N+1) V\k a ,n a ) . (16) 
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We will isolate the dependence of £ on the density fluctuations and determine the extent of their impact on the 
suppression of the fringe contrast. This can be accomplished analytically in the regime in which the energy of the 
probe, Eq = ft 2 ko / (2m) , is large compared to the bandwidth of the optical lattice, but not sufficiently large to excite 
atoms out of the lowest band, so that we continue to have single band dynamics. The optical lattice potential along 
each spatial direction is Vl (x) = Vl sin 2 (x) , where the depth of the lattice is conveniently specified in units of the 
photon recoil energy E r — h 2 k L / '(2tot). fez, is the laser wave number, and mr is the mass of an atom in the target. 
For a typical lattice depth of Vl = 15E r , the width of the lowest band, determined by the energy difference between 
the center and the edge of the first Brillouin zone, is found from the exact solution for the optical lattice eigenvalue 
equation to be 0.03E r . Likewise, the gap between the first and second bands, determined by the energy difference 
between these bands at the edge of the first Brillouin zone, is 6.28E r . It is, therefore, readily possible to employ a 
probe with energy much larger than the lowest band's width, which is nonetheless insufficient to bridge the band gap. 
In this situation, we may neglect the loss of energy by the probe. The second order term in the expansion of the 
diagonal T-matrix element that appears in Eq. (10) is 

(ho, n | VG { N+1) V \ko, n ) = ^ J d 3 r dV j*o-(r>-r) ^ - (r) q(n+i) ^ r ,. Eq + ~ (p/) ^ _ 

The matrix element in this expression has a clear physical interpretation. The probe wave function is scattered first 
at r', with a strength proportional to the density of atoms in the target at this location. The scattered probe then 
evolves freely according to the Green's function from r' to r, where a second scattering event takes place. As we are 
computing the coherently scattered projectile wave function, we only include the contribution to this process which 
leaves the target untouched, in the state \n ). We then integrate over all such two-point scattering paths. 

Although the target must be left untouched, there may be transient excitations between the two scattering events. 
Formally, this can be seen by inserting a complete set of target states, ^ n \n) (n\, into the transition matrix element 
in Eq. (17). Note that the eigenstates of the target, |n), are also eigenstates of the many body Green's function, 

J2G N+1) (r,r' ]Eo + E no ) \n) (n| = £ <r| — — 1 \r') \n) (n\ . (18) 

n n E Q + (E na - E n ) - Hp +l£ 

The coefficients in this expansion are the usual, single particle free Green's functions at the shifted energy s n = 
Eo + (E m - E n ), 

G (r, r'; e n ) = (r\ 1 \r') . (19) 

e n - Hp + it 

Physically, each term corresponds to an intermediate inelastic scattering channel, in which the projectile has trans- 
ferred energy E n — E no to the target. The largest possible energy transfer due to scattering a single atom in the target 
within the lowest band is the bandwidth. As we are considering the situation in which the initial probe energy is large 
compared to the bandwidth of the lowest band, the transient shift in its energy will be comparatively small, so that 
e n w Eo. The single particle Green's function that appears in the expansion in Eq. (18) is, therefore, approximately 
independent of the sum over target states. Inserting this expansion into the second-order correction to the T-matrix 
element given in Eq. (17) yields 



(ko, n | VG ( Q N+1) V |fco, no) » ^ / d 3 r d 3 r' e ik °< r '-^G Q (r, r'; E Q ) (n | h(r) h(r') \n ) . 



(20) 



In arriving at this expression, we have made use of the fact that the projectile loses very little energy as it passes 
through the target, and that interactions with individual atoms in the target are predominantly s-wave scattering, due 
to the slow relative velocities. The amplitude for each individual double scattering sequence appropriately depends 
on the density-density correlator between the two scattering points. 

As we seek the relationship between the on-site particle number fluctuations in the lattice and the decay of the 
contrast of the interference fringes, we must express the density-density correlator in Eq. (20) in terms of the on-site 
field operators. The field operator for the atoms in the target can be expanded in terms of the Wannier function for 
the lowest band of the lattice, w(r), as 

N L 

^(r) = Y,a j w{r-R j ), (21) 
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where the operator dj annihilates a particle at the j th lattice site. The density n(r), in terms of this expansion, is 
then given by 

h(r)= ^2a]a k w*(r-R j )w(r-R k ). (22) 
j,k=i 

For the lattice depths that we are considering, Vl ~ 15E r , we may employ the tight-binding approximation and 
neglect the off-diagonal terms due to the small overlap between Wannier functions centered on different lattice sites. 
The relationship between the density at a position r and the density at a site, j, is then approximately 

N L 

nW^^-Kr-^l 2 . (23) 

Substituting this into the density-density correlator will allow us to explicitly extract the dependence on the fluctua- 
tions, 

N L 

(no| n(r) n(r') |n ) « ^ (n,n fe )|w;(r - fl,)| 2 |tu(»-' - #fc)| 2 - (24) 
j,k=i 

When the number of atoms N in the target is much larger than f , a condensate of N atoms in the ground state of the 
lattice is well approximated by a coherent state with an average of N atoms. In the weakly interacting case, therefore, 
the target ground state is appropriately described as a product of coherent states at each lattice site with an average 
of n atoms per site [3, p. 906]. In the strongly interacting case, the target is a product of Fock states with n atoms per 
site. In either case, we may simplify the cross correlators between different lattice sites, {hjhk) — {hj){hk) (j ^ k), so 
that these correlators are independent of the many body phase. For a uniform lattice, the average density does not 
vary from site to site, so that (nj) — n — N/Nl. The diagonal terms in the sum yield a dependence on the second 
moment of the on-site number distribution, which is also independent of the particular index, j. The second moment 
can be re-expressed in terms of the on-site fluctuations, u, as a 2 — (n 2 ) — n 2 . The contribution of the square of the 
average density may be combined with the off-diagonal terms to give a sum over all pairs of sites, leaving 

N L N L 

(n a \h(r)h(r')\n Q ) ^ a 2 Y,Hr ~ Rj)\ 2 \w(r' ~ Rj)\ 2 +n 2 £ \w(r - Rj)\ 2 \w(r' - R k )\ 2 . (25) 
j=i j,k=i 

The second term in the density-density correlator, proportional only to the average density, will not vary as the 
interaction strength is increased across the transition between the superfluid and Mott insulating phases; however, 
the first term, proportional to the density fluctuations, will vanish as the interaction strength between the atoms in 
the lattice increases, becoming zero in the Mott insulator. The appearance of fluctuations in the diagonal T-matrix 
element will lead to an additional loss of contrast in the interference fringes. Inserting this expression for the density- 
density correlator into the expression for the second order contribution to the T-matrix element given in Eq. (20) 
yields 

(k ,n \VG { N+1) V\k ,n ) = ^ffi F< 2 > (fc ) + A, (26) 

where A refers to the contribution of the term proportional to n 2 . This quantity is independent of the interaction 
strength as the average density of atoms per lattice site is fixed for any U/J and will therefore not contribute to 
the variation of the interference fringe contrast. The function F^ 2 \k ) gives the amplitude for all possible two-point 
scattering events, 

F {2) {k ) = j d 3 rd 3 r' e ik °< r '-^ \w{r)\ 2 G (r, r'; E ) \w(r')f . (27) 

The imaginary part of F^(ko) is calculated in Appendix A for the concrete situation in which the lattice is one- 
dimensional and oriented perpendicular to the incident probe beam (see Fig. 1). It is then given by 
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where s — Vl/E t is the depth of the lattice in units of the photon recoil energy. Using this, we can explicitly 
write down the factor that modifies the contrast of the interference fringes due to the quantum mechanical density 
fluctuations in the lattice, 

_ / pal a 2 sV* vt f ko/kL \\ f0Q , 

Here we note that the functional dependence of the contrast on the wave number of the probe and the depth of the 
lattice is in terms of the dimensionless quantity K — (fco//cL)/(v / 2s 1 / 4 ). Employing K, the factor giving the reduction 
of the contrast due to the fluctuations in the lattice takes the form 

Cocexpf-^=^-^Erf(if)V (30) 



40F n K 

For a typical lattice depth, s — 15, and for the range of probe wave numbers which avoid interband excitations at 
this lattice depth, ko/k^ < V6.28, the variation of the exponent with K is minor, and the function, K~ l Eii(K) is of 
order 1. The remaining dependence of the contrast is such that it decays exponentially with increasing fluctuations 
in the on-site atom number and with increasing column density of the sample. 

The attenuation of contrast given in Eq. (30) can be directly obtained from experimental measurement of the 
interference fringes as the ratio of the contrast of the fringes at arbitrary values of the parameter U/ J to the contrast 
deep in the Mott insulator regime, where a 2 — > 0. Although it is common to alter the interaction strength U/J by 
adjusting the depth of the lattice, here we consider adjusting U/J by manipulating the scattering length of the lattice 
atom collisions through a Feshbach resonance, so that we retain a constant lattice depth, and K remains constant. 

We have shown that the quantum mechanical density fluctuations at individual lattice sites have a decohering effect 
on a matter wave passing through a Mach-Zender interferometer, and that the contrast in the interference fringes 
is exponentially suppressed by increasing fluctuations. It was necessary to employ an interferometric arrangement 
to explicitly extract the coherent component of the scattering, and in this way, we have definitively verified the 
enhancement of incoherent scattering due to these density fluctuations. In the following section, we will exploit 
these results to identify the influence of the density fluctuations on the scattering cross section and demonstrate the 
conditions under which the fluctuations are directly observable therein. 



III. DENSITY FLUCTUATIONS AND INELASTIC SCATTERING 



We begin our analysis of the differential scattering cross section with a general result for scattering from a many 
body system in the first Born approximation. The target is initially in the ground state, |n ). The cross section is 
then given by [10, 22] 

2 

(31) 

where n = k — k is the difference between the incident probe wave vector k and the outgoing probe wave vector k. 
The magnitude of k is potentially reduced from the incident probe wave number, due to a transfer of energy from the 
probe to the lattice, and is determined by H 2 k 2 /(2m) — h 2 k 2 /(2m) + (E no — E n ). The outgoing probe wave vector 
points in the direction 9 of the detector (see Fig. 1). 

Bragg peaks, which correspond to a coherent sum of waves scattered from individual atoms in the lattice, appear 
in the elastic scattering cross section [10]. Additionally, we showed in Sec. II that increased density fluctuations 
correspond to increased incoherent scattering of the probe. This suggests that the on-site density fluctuations will 
influence the scattering cross section through inelastic channels, in which the probe delivers energy to the atoms in 
the lattice. This intuition motivates us to examine the cross section under the same energetic conditions as in our 
analysis of the interference pattern in Sec. II, described immediately following Eq. (16). There, the contribution of 
single band inelastic channels was emphasized by choosing an initial probe energy sufficient to allow excitation to all 
of the modes of the lowest band of the lattice, but which did not exceed the splitting between the lowest band and 
the first excited band. 

The transition amplitude between the ground state, no, and an excited state, n, due to a momentum boost, which 
is given by the matrix element in the expression for the cross section in Eq. (31), is non-zero in the first Born 
approximation only for final target states in which at most a single atom in the target is displaced out of the ground 
state of the lattice. The energy of the final target modes that we must consider, E n , differ, therefore, from E by the 



da 
dtt 



E - E 
fi 2 fc2/(2m) 



N 

HE' 

3 = 1 



K) 
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energy of a single particle excitation. Thus, the conditions are viable to enforce the condition E n — E na <C fi 2 ^/(2m) 
for all accessible final target states. At such high probe energies, the factor under the square root in Eq. (31), which 
weights the contribution of the scattering into the final target mode, n, is nearly uniform and approximately unity. 
This factor appears also in k; however, in the high-energy approximation, we take k to be independent of the energy 
transfer. This is equivalent to the so-called static approximation [22]. With these simplifications, the cross section 
becomes 



da 



N 



3=1 



5>) (n 




K) . 



(32) 



In the subspace of the lowest band of the lattice, the sum over projections onto individual target modes, n, is an 
identity and can be removed. The cross section can then be expressed in terms of the density-density correlator, given 
in Eq. (25), by rewriting it in second quantized notation. This is accomplished via the identification of J^fcLi e tK, ' rk 
with Jd 3 r e lK ' r h(r), which leads to the expression for the cross section, 



da 



: J d 3 rd 3 r' e -^-{r-r>) a( P )ft( P ') \ Ho ) . 



(33) 



We may expand the density n(r) according to Eq. (23), in terms of on-site densities and the Wannier function of 
the lowest band. When the number of lattice sites is sufficiently large that the effect of the edges may be neglected, 
the cross section separates into two factors. One of which is a smooth background determined by the shape of the 
Wannier function, which reflects the depth of the lattice. The other depends only on the distribution of atoms in the 
lattice. Thus for Nl ~> 1, we find 



a 



1 da 
~2~dtt 




in(R h -Rj) 



1 



\w(r)f 



(34) 



A completely filled lattice, i.e., N > Nl, with a large number of lattice sites will have a correspondingly large number 
of atoms, so that the superfluid ground state is well-described by a coherent state with an average of N atoms. Then 



we have as before that (hjfik) 



(j ^ fc) and {hfj = a 



The expression for the cross section, valid in both 



the superfluid and Mott insulating phases, is given by, 



1 da 
aldtt 



e -iK-R, 

3=1 



\ 



+ a 2 N L 



J 



I 



d 3 r e- iK - r \w(r){ 



(35) 



In this expression, the first term gives elastic Bragg peaks that scale as the square of the number of atoms in the 
lattice. The remaining term gives the inelastic cross section, which contributes a smooth background. In the high 
probe energy regime, the inelastic background is directly proportional to the fluctuations. Its amplitude will therefore 
be reduced as the interaction strength between the atoms in the lattice is increased relative to the tunneling matrix 
element. The fluctuations are shown as a function of the parameter U / J in Fig. 3 for a lattice with N = 6 atoms and 
Nl = 6 lattice sites. This is sufficient due to the rapid convergence of the shape of the decay of the fluctuations with 
increasing lattice size [24]. The sum over exponentials in the first term yields Bragg peaks located at positions for 
which k is a reciprocal lattice vector. These arise due to the coherent overlap of the scattered waves from individual 
atoms in the lattice. Our result for the cross section recreates the expected behavior, in which, deep in the Mott 
insulator phase, with a 2 — > 0, no inelastic background remains, whereas the superfluid phase, in which a 2 w n, yields 
an inelastic background proportional to the number of atoms in the lattice, N (see Fig. 4). 

The final factor, which is the Fourier transform of the Wannier function, leads to an overall envelope with an 
approximately Gaussian shape. In the region between the extremes of very weakly or very strongly interacting 
atoms in the lattice, the fluctuations can be directly observed by measuring the scale of the inelastic background. In 
particular, the difference between the cross section for an arbitrary value of U/J and the cross section deep in the 
Mott insulator regime will eliminate the elastic Bragg peaks, and leave only a Gaussian background that scales as 
o 2 Nl- This result is intuitively satisfying in that we can trace the impact of the density fluctuations to the incoherent 
part of the scattering. This is made explicit through our analysis of the interferometer, and likewise is consistent with 
our expectations for the scattering cross section, in which we associate the inelastic part of the cross section with 
incoherent scattering processes. 
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FIG. 3. The on-site number fluctuations, a, plotted as a function of the parameter U/J, for a lattice with N = 6 atoms and 
Nl — 6 lattice sites and a depth of s = 15 times the photon recoil energy. 
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FIG. 4. The cross section given by Eq. (35) for scattering a matter wave from the atoms in a one-dimensional lattice, oriented 
perpendicular to the incident probe wave vector, as a function of the angle 9 from the forward direction. The first term in the 
cross section, which gives Bragg peaks, are computed for 100 atoms in 100 lattice sites. The second term, giving the inelastic 
background, is plotted using the values of a given in Fig. 3 and a tunneling matrix element J = 0.0047_E r - The plot in light 
gray corresponds to U/J = 0, in which the target is entirely superfluid. The small deviation from N is a result of using the 
exact, non-interacting ground state, rather than the coherent state approximation. The cross section for scattering from a 
target with U/J = 20 is shown in dark gray. 



IV. CONCLUSIONS 

We have shown that matter wave scattering is susceptible to the purely quantum mechanical fluctuations in on-site 
atom number in an optical lattice. The fluctuations lead to incoherent scattering of the probe atom, which can be 
quantified and observed by the loss of contrast of the interference fringes in a Mach-Zcndcr interferometer. The 
fringes form due to the interference between the scattered matter wave and the unscattered, reference arm of the 
interferometer. To first order in the interaction strength, Vb, the interference pattern is shifted due to a quantum 
mechanical phase proportional to the density of atoms in the lattice. Decoherence occurs in the interferometer as 
a second order effect. The second order term in a Born series expansion of the many body scattering T- matrix, 
proportional to Vq, corresponds physically to double scattering of the probe within the cold atom sample. When this 
multiple scattering is taken into account, the scattered matter wave becomes sensitive to the local density fluctuations 
in the target, which cause decoherence of the probe atom. 

In the regime in which the scattered probe atom has an energy that is large compared to the bandwidth of the lowest 
band of the optical lattice, the fluctuations can be directly obtained from the amplitude of the inelastic scattering 
background. Whereas the elastic part of the scattering cross section corresponds to coherent scattering, the inelastic 
part contains the incoherent contribution, in which atoms in the lattice are excited out of the ground state. This 
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incoherent process exhibits a simple and direct dependence on the density fluctuations, providing a method to directly 
observe them. 
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Appendix A: Derivation of F* 2 '(fc ) 

We will evaluate the two-point scattering function given in Eq. (27), 

(fc ) = J d 3 r d 3 r' e ik °< r '-^ \w{r)\ 2 G (r, r'; E ) \w(r')f , (Al) 

explicitly for the case of a one-dimensional lattice, in which the atoms are tightly confined in the y and ^-directions, 
and the wave vector ko of the incident probe atom points along the z-direction. Strictly speaking, the range of the 
integration is over the length of the lattice; however, the locality of the Wannier function allows us to extend the 
integration to all space. For a harmonic confinement that is deep in the y and ^-directions, the Wannier function 
approaches |w(r)| 2 « |w(a;)| 2 5{y) S(z). The one-dimensional case, with the substitution u = x — x' , therefore yields, 

r+cc ik \u\ r+oo 

i " (2) (feo) = -^ TI / du j : / dx'\w{u + x')\ 2 \w(x')\ 2 . (A2) 
2tt^ J_ oa \u\ J_ oc 

We can use the harmonic approximation for the Wannier function in a lattice of depth s = Vo/E r , 

w(x) « (jljl.ErVij ' exp (-^E r ^x 2 /2^ . (A3) 
This allows us to compute the integral over x' in Eq. (A2) analytically. The imaginary part of F( 2 )(fc ) is then given 

by 

Im(F( 2 )(fc )) = - m "iX^ 2 S -^e~^^ 2 . (A4) 

Notice that we have taken advantage of the evenness of the integrand to reduce the range of the integration to [0, oo) 
and thereby eliminated the need for the absolute value operation on u. The remaining integral over u is then given 

by 
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